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Abstract 

Existence of a solution to the quasi- variational inequality problem 
'^ . arising in a model for sand surface evolution has been an open problem 

for a long time. Another long-standing open problem concerns deter- 
mining the dual variable, the flux of sand pouring down the evolving 
sand surface, which is also of practical interest in a variety of appli- 
cations of this model. Previously, these problems were solved for the 
special case in which the inequality is simply variational. Here, we 
introduce a regularized mixed formulation involving both the primal 
(sand surface) and dual (sand flux) variables. We derive, analyse and 
compare two methods for the approximation, and numerical solution, 
of this mixed problem. We prove subsequence convergence of both 
approximations, as the mesh discretization parameters tend to zero; 
f— ^ I and hence prove existence of a solution to this mixed model and the 

CN ' associated regularized quasi-variational inequality problem. One of 

these numerical approximations, in which the flux is approximated 
by the divergence-conforming lowest order Raviart-Thomas element, 
^ ' leads to an efficient algorithm to compute not only the evolving pile 

^ ■ surface, but also the flux of pouring sand. Results of our numerical 

experiments confirm the validity of the regularization employed. 
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1 Introduction 

Let a cohesionless granular material (sand), characterized by its angle of 
repose a, be poured out onto a rigid surface y = wo{x), where y is vertical, 
xeQg M.'^, d = 1 ot 2, and 1] is a domain with boundary dQ. The support 
surface wq e Wq''^{Q) and the nonnegative density of the distributed source 
/ G L^(0,T; L'^{fi)) are given. We consider the growing sandpile y = w{x,t) 
and set an open boundary condition w\gn = 0. Denoting by q{x,t) the 
horizontal projection of the flux of material pouring down the evolving pile 
surface, we can write the mass balance equation 

The quasi-stationary model of sand surface evolution, see Prigozhin fi6\ \T7\ 
HE], assumes the flow of sand is confined to a thin surface layer and directed 
towards the steepest descent of the pile surface. Wherever the support surface 
is covered by sand, the pile slope should not exceed the critical value; that is, 
w > Wq ^ \'V_w\ < ko, where k^ = tana is the internal friction coefficient. 
Of course, the uncovered parts of support can be steeper. This model does 
not allow for any flow on the subcritical parts of the pile surface; that is, 
|Vw| < ko =^ 1 = O- These constitutive relations can be conveniently 
reformulated for a.e. (x, t) E fl x (0, T) as 

\Vw\<M{w) and M{w) \q\ +^10 .q = 0, (1.2) 

where for any xEQ 

M(w)(x) := I ^' «^(^) > «^o(x), .^3. 

I va.ax{ko,\V_Wo{x)\) i^ix) < Wo{x). 

Let us define, for any r] G C{Q), the closed convex non-empty set 

K{7]) := {if G Wo'°°(^) : iVy^l < M{r]) a.e. in n} . (1.4) 

Since M{w) \q\ + Vip . g > for any ip G K{w), we have, on noting (II. 2p . that 
w G K{w) and V((/3 — w).q>0. A weak form of the latter inequality is: for 
a.a. t G (0,T) 

V.q{w-if)dx>0 \f(peK{w). (1.5) 

Jn 



Combining f ll.Sp and fll.ip yields an evolutionary quasi- variational inequality 
for the evolving pile surface: Find w G K{w) such that for a.a. t G (0,T) 

(dw \ 
-7-- f](v-w)dx>0 W^eKiw). (1.6) 

dt J 

Assuming there is no sand on the support initially, we set 

w{;0)=Wo{-). (1.7) 

A solution w to this quasi- variational inequality problem, (II. 6p and (II. 7p . 
if it exists, should be a monotonically non- decreasing function in time for 
any / > 0, see §3 in Prigozhin [TH]. However, existence and uniqueness of 
a solution has only been proved for support surfaces with no steep slopes; 
that is, |Vwo| < ko, see Prigozhin [161 IIH]- In this case K{w) = K : = 
[(f G Wq'°°{Q) : \Vsp\ < ko a.e. in fi^ and the quasi- variational inequality 
becomes simply a variational inequality. Independently, the variational in- 
equality for supports without steep slopes has been derived and studied in 
Aronson, Evans and Wu [IJ as the p — )■ oo limit of the evolutionary p- 
Laplacian equation. 

The quasi-variational inequality (II. 6p can, of course, be considered not 
only with the initial condition (II. 7p . However, if w(-, 0) = Wo{-) > Wo{-) and 
Wq does not belong to the admissible set K{wq), an instantaneous solution 
reconstruction takes place. Such discontinuous solutions, interpreted as sim- 
plified descriptions of collapsing piles with overcritical slopes, were studied 
in the variational inequality case in Evans, Feldman and Gariepy [10], and 
Dumont and Igbida [8J. Since we assumed the initial condition (II. 7p and, 
obviously, Wq G K{wo), one could expect a solution continuously evolving in 
time. However, for the quasi-variational inequality with the open boundary 
condition w\gQ = 0, an uncontrollable influx of material from outside can 
occur through the parts of the boundary where V_Wq . z/ > k^, where u_ is the 
outward unit normal to dfl. This makes the solution non-unique and, pos- 
sibly, discontinuous. Such an influx is prevented in our model by assuming 
that 

Vwo ■i!.< ko on do,, (1.8) 

which implies that V_w .iJ_<kQon dQ for t > 0. 

For the variational inequality version of the sand model, equivalent dual 
and mixed variational formulations have recently been proposed; see, e.g.. 



Barrett and Prigozhin |1] and Dumont and Igbida [7]. Such formulations 

are often advantageous, because they allow one to determine not only the 
evolving sand surface w but also the surface flux q, which is of interest too 
in various applications; see Prigozhin [l5l[T7], and Barrett and Prigozhin |6]. 
In such formulations, and this is their additional advantage, the difficult to 
deal with gradient constraint is replaced by a simpler, although non-smooth, 
nonlinearity. 

Here we will also use a mixed variational formulation of a regularized 
version of the growing sandpile model involving both variables. Instead of 
excluding the surface flux q from the model formulation, as in the transition 
to fll.6p above, we now note that the first condition in fll.2p holds if for a.e. 
(x,t) G fix (0,T) 

M{w)\v\+Yw.v>0 (1.9) 

for any test fiux v_. Hence we can reformulate the conditions (11. 2p for a. a. 
te (0,T) as 

/ [M{w){\v\-\q\)-w'V.{v-q)]dx>0 (1.10) 

Jn 

for any test fiux v, and consider a mixed formulation of the sand model as 
( ITT]) and (ILTOJ) . 

The quasi- variational inequality (I1.6P is a difficult problem; in particular, 
due to the discontinuity of the nonlinear operator M, which determines the 
gradient constraint in (II. 4p . Furthermore, the natural function space for 
the fiux q is the space of vector- valued bounded Radon measures having L^ 
divergence. If q is such a measure, the discontinuity of M{w) also makes it 
difficult to give a sense to the term J^ M{w) \q\ dx in the inequality (ll.lOp of 
the mixed formulation. 

In this work we consider a regularized version of the growing sandpile 
model with a continuous operator M^ : C{Q) — > C{Q), determined as follows. 
For a fixed small e > 0, we approximate the initial data Wq G Wq'°°{Q) by 
Wq G WQ''^{Q)r\C^{Q), and M(-) by the continuous function M^{-) such that 
for any x G fi 

M,{r)){x) 



r]{x) > Wq{x) + e, 
r]{x) E [wq{x] 

ki{x) '■= max(A;o, |Vti'Q(a;)|) ri{x) < Wq{x) 



klix) + {ko - k\{x)) ( "^^^ ^"^"^-^ ) r]{x) E Hi^),w'o{x) + e] 



1.11) 



We note that M^ is such that for all r/i, 772 ^ C(i7) 



|M,(r7i) - M,(r/2)|o,oo,Q < -^^ \vi - mWo^n , (1.12) 



where 



kl^:= max kl{x). (1.13) 

xen 

In addition, it follows for any xeVI that 

r)i{x) > r)2{x) => 0<ko< M,{r)i{x)) < M,{r,2{x)) < fc^(x) . (1.14) 

We note that the analysis of the sand quasi- variational inequality problem 
studied in this paper is far more involved than that of the superconductivity 
quasi- variational inequality problem studied by the present authors in [6]. 
In the superconductivity context, M : R -> [Mq, Mi] C M with Mq > 0. 
In |6], we exploit the fact that |Vu;(a;.)| < M{w{x)) can be rewritten as 
\Y[Fiw{x))]\ < 1 for all x E n, where F'(-) = [M(-)]-i and F(0) = 0. 
Clearly, such a reformulation is not applicable to M(-), (11. 3p . or M^{-), (11. lip . 

In addition, we note that in the very recent preprint by Rodrigues and 
Santos [in] an existence result can be deduced for the primal quasi- variational 
inequality problem (ll.6|) for a continuous and positive M(-), such as M^[-), 
and / E W^'°°{n X (0,r)). Assuming w° G K{w°) n Co(n), they show that 
w E L~(0,T; Wo'°°(fi)) H W^'°°{0,T; [Co{n)Y). Their proof is based on the 
method of vanishing viscosity and constraint penalization. 

The outline of this paper is as follows. In the next section we introduce 
two fully practical finite element approximations, (Q^'"^) and (Q^'"^), to the 
regularized mixed formulation (II. ip and (ll.lOp . where M(-) is replaced by 
Mir{-), and prove well-posedness and stability bounds. Here h and r are the 
spatial and temporal discretization parameters, respectively. In addition, 
r > 1 is a regularization parameter in replacing the non-differentiable non- 
linearity I ■ I by the strictly convex function - | ■ |^. The approximation (Q^'^) 



is based on a continuous piecewise linear approximation for w and a piece- 
wise constant approximation for g, whereas (Q^'^) is based on a piecewise 
constant approximation for w and the lowest order Raviart-Thomas element 
for q. In Section |3] we prove subsequence convergence of both approximations 
to a solution of a weak formulation of the regularized mixed problem. This is 
achieved by passing to the limit h —> Q first, then r — )■ 1 in the case of (Q^'"^), 
and finally r — )■ 0. In Section HI we introduce iterative algorithms for solving 
the resulting nonlinear algebraic equations arising from both approximations 
at each time level. Finally, in Section [5] we present various numerical exper- 
iments. These show that only the approximation (Q^'^) leads to an efficient 
algorithm to approximate both the surface w and the flux q. 

We end this section with a few remarks about the notation employed 
in this paper. Above and throughout we adopt the standard notation for 
Sobolev spaces on a bounded domain D with a Lipschitz boundary, denoting 
the norm of W^'^{D) (^ G N, s G [l,C)o]) by ||.||^,s.z) and the semi-norm by 
I ■ V,s,D- Of course, we have that | ■ |o,s,d = || ■ ||o,s,d- We extend these norms 
and semi-norms in the natural way to the corresponding spaces of vector 
functions. For s = 2, W^''^{D) will be denoted by H^[D) with the associated 
norm and semi-norm written as, respectively, || ■ H^^/) and | ■ l^^^). We set 
Wl'\D) ■= {r] G W^''{D) ■ T] = OondD}, and H^{D) = W^'^D). We 
recall the Poincare inequality for any s G [1, oo] 

\v\o,s,D < C,{D) \Vr]\o,s,D yv^ W^'^iD) , (1.15) 

where the constant C^,{D) depends on D, but is independent of s; see e.g. page 
164 in Gilbarg and Trudinger [13]. In addition, \D\ will denote the measure 
of D and (■, ■)d the standard inner product on L'^{D). When D = Q, for 
ease of notation we write (■, ■) for {■,-)n- 

For m G N, let (i) C"^{D) denote the space of continuous functions with 
all derivatives up to order m continuous on D, (ii) C^{D) denote the space 
of continuous functions with compact support in D with all derivatives up to 
order m continuous on D and (iii) C^{D) denote those functions in C"^{D) 
which vanish on dD. In the case m = 0, we drop the superscript for all 
three spaces. 

As one can identify L^{D) as a closed subspace of the Banach space of 
bounded Radon measures, Ai{D) = [C{D)]*, i.e. the dual of C{D); it is 



convenient to adopt the notation 

/_l/^l = \\f^\\M(D) ■= sup {fi,v)c(D) < OO, (1.16) 

l»?lo,oo,D<l 

where (■, ■)c(d) denotes the duahty pairing on [C{D)Y x C{D). 
We introduce also the Banach spaces for a given s G [1, cxd] 

V\D) := {v e [L'{D)]^ -.V-ve L\D)} , (1.17a) 

V^{D) :={ve [M(D)Y:V.veL\D)}. (1.17b) 

The condition V_.vE L'^{D) in (I1.17a[ b) means that there exists u G L'^{D) 
such that {v,'V(J))q(jj\ = —{u,(f))D for any G Cq{D). 

We note that if {/i„}„>o is a bounded sequence in Ai{D), then there exist 
a subsequence {finj}nj>o and a /i G A^(i5) such that as rij — )■ oo 

fin, -> yu vaguely in M{D); i.e. (/i„^, - /i, ?7)c(d) ^0 '^ ^7 ^ C(L') . 

(1.18) 

In addition, we have that 

liminf / |/i„.| > / |/i| ; (1-19) 
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see e.g. page 223 in Folland [T2]. 

We note that V-^{D) and Vf{D) for s _G [1,2) are not of local type; 
that is, V G V-^{D) [Vf{D)] and (p G C°°(U) does not imply that (pv e 
y_^{D) [Vf{D)], see e.g. page 22 in Temam [22]. Therefore, one has to avoid 
cut-off functions in proving any required density results. If Q is convex, which 
we shall assume for the analysis in this paper, then it is strictly star-shaped 
and one can show, using the standard techniques of change of variable and 
mollification, that 

[C°°(n)]2 is dense in Vf{n) if sG(1,oo). (1.20) 

Moreover, for any v G V/^{il), there exist {Vj}j>i G [C°°(f2)]'^ such that 

Vj — )■ V vaguely in [M{Q)]''' as j — )■ oo , (1.21a) 

Y . t; J -)■ V . t; weakly in L^(fi) as j -)■ cx) , (1.21b) 



limsup / p\Vj\dx= / p\v\ (1.21c) 

for any positive p G C(r2). We briefly outline the proofs of (]1.21al -c). With- 
out loss of generality, one can assume that Q is strictly star-shaped with 
respect to the origin. Then for v defined on Q and ^ > 1, we have that 
^el^) — ^(^~^2^) is defined on Qq := 6 fl D fl. Applying standard Friedrich's 
moUifiers J^ to Vg, and a diagonal subsequence argument yield, for ^ ^ 1 
and e — )■ as j — )■ oo, the desired sequences {v}j>i demonstrating ( ]1.20p if 
v_ G Vf{^) and satisfying f ll.21af c) iiV_^{il); see e.g. Lemma 2.4 in Barrett 
and Prigozhin |5], where such techniques are used to prove similar density 
results. 

We recall the following Sobolev interpolation theorem, see Theorem 5.8 
in Adam and Fournier [2]. If 77 G W^'^{D), with s > d, then 77 G C{Q) with 
the embedding being compact; and moreover, 

ds 

\v\o,oo,D<C{s,D)\\ri\\l^^o\r]\l-j^ with a = ^^ ^ ^, _ ^. G (0, 1) . 

(1.22) 

We recall also the Aubin-Lions-Simon compactness theorem, see Corollary 
4 in Simon [2^. Let Bq, B and Bi be Banach spaces, Bi, i = 0,1, reflexive, 
with a compact embedding Bq ^-^ B and a continuous embedding B ^-> Bi. 
Then, for a > 1, the embedding 

{r/GL-(0,T;i3o):^GL"(0,T;i3i)}-^C([0,r];i3) (1.23) 

is compact. 

Finally, throughout C denotes a generic positive constant independent of 
the regularization parameter, r G (1, 00), the mesh parameter h and the time 
step parameter r. Whereas, C{s) denotes a positive constant dependent on 
the parameter s. 



2 Finite Element Approximation 

We make the following assumptions on the data. 

(Al) Q C W^, d = 1 OT 2, is convex having a Lipschitz boundary dQ with 
outward unit normal v_. f & L'^{0,T]L'^{VL)) is a nonnegative source, and 



M^(-) is given by (ITTT]) . In addition, the initial data w^ G Wo''^{n) n C^((]) 
is such that V_Wq .iy_< k^. 

For ease of exposition, we shall assume that fi is a polygonal domain to 
avoid perturbation of domain errors in the finite element approximation. We 
make the following standard assumption on the partitioning. 

(A2) Vt is polygonal. Let {T'^}h>o be a regular family of partitionings of VL 
into disjoint open simplices a with h^j := diam(cr) and h := maXg-gT-h h^j, so 
that VL = Ug-gT-ha. 

Let Uq^ be the outward unit normal to da, the boundary of a. We then 
introduce the following finite element spaces 

S'':={r]''eL^{n):r]''l=a,eM WaeT^}, (2.1a) 

5^0 := {^' e L°°(fi) : r^" U= a^ G M>o '^aeT'^}, (2.1b) 

S^:={T/'e[L^in)f:rf\^=a^eR'' \/ a E T'' } , (2.1c) 

U'' := {r/^ G C(n) : r/^ U= o^ . X + 6,, o^ G M^ 6, G M V a G T^} , 

(2.1d) 

[/q^ := ^7'^ n Vro'"(^) , (2-le) 

Y!" := {t;^ G [L^{n)Y : t;'' U= a, + b^x, g^ eR'^, b, eR V a G T^' 

and (/ U -^'^ U') -^a = on (9cT n da' \/ a, a' e T^} . (2.1f) 

Here V_ is the lowest order Raviart-Thomas finite element space. 

Let 71^ : C{Q) — )■ U^ denote the interpolation operator such that 
Ti^r]{x,j) = ri{Xj), j = 1 -^ J, where {xj}^^^ are the vertices of the parti- 
tioning T^. We note for 171 = and 1 that 

1(1 - 7r'')vU,s,a < C h^-"" \7]\2,s,. yaer\ for any sG[l,oo], (2.2a) 
lim 11(7 - 7r'')r]\\m,oo,n = y 7] e C^m ; (2.2b) 

where I is the identity operator. Let P^ : [L^{ri)Y — > S_^ be such that 

P!'v\a= 7^ [ vdx yaeT^. (2.3) 

We note that 

\P^v\o,s,a < \v\o,s,a ^ H ^ [L' (a)]'' , sG[1,oo], WaeT^ (2.4a) 



lim I \v\ - \P%\ |o,oo,n < lim \v - P^|o,oo,c = V t; G [C{n)Y . (2.4b) 

Similarly, we define P^ : L^{Q) -^ S^ with the equivalent to fl2.4a[ b) holding. 
In addition, we introduce the generalised interpolation operator J^ : 
[Vr^'*(n)]'^ -)■ y_'', where s > 1, satisfying 

v-L'\)-Es,ads = i = 1^3, VaeT''; (2.5) 

di(T 

where da = Uf^idicr and z/g.^ are the corresponding outward unit normals 
on dia. It follows that 

{V.{v-tv),v'') = ^v'^^S''. (2.6) 

Moreover, we have for all a G T'^ and any s & (1, C)o] that 

Will - \L^ll\\o,s,a < \V~ L^v\o,s,a < C h„\v\i^s,a and \l^v\l,s,a < C \v\i^s,a , 

(2.7) 

e.g. see Lemma 3.1 in Farhloul [11] and the proof given there for s > 2 is 
also vahd for any s G (1, oo]. 

We introduce (?7,x)^ := Eaer'^(^'^)^ ^^^ 

d+l 



(j,-r_L « 

^1^1 $^r/(a;J)x(a;J)= Avx]dx 
i=i -^'^ 



Vr/, xGC(a), V^gT"; (2.8) 

where {xj} jt} are the vertices of a. Hence (r], x)'^ averages the integrand r] x 
over each simplex a at its vertices, and is exact if 77 x is piecewise linear over 
the partitioning T'^. We recall the well-known results that 

\vXa < IvX ■■= {V\ rl^f < (rf + 2) \riXn ^ "<t ^U\ (2.9a) 



|(^',x')-(r]',x')1 = |((/-7r^)(r/'^x'),l)| < |(/-7r^)(r/'^x')|o. 



i,n 



<C/i|r/^|o,n|x'|i,n Vr/^x^ef/^ (2.9b) 
where we have noted 02.2ap . 
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In order to prove existence of solutions to approximations of fll.lOl) . we 
regularise the non- different iable nonlinearity | ■ | by the strictly convex func- 



tion 



for r > 1. We note for all a, b eM. that 



1 dial'' 
r daj 



I ir-2 

\a\ aj 



ar a.{a-b) > i [ |ar - |6|' 



(2.10) 



Similarly to fl2.9ap . we have from the equivalence of norms and the convexity 
of I ■ 1^ for any r > 1 and for any v^ eY_ that 

C(|z;'^r,l)^< /"|t;^rdx < (|/r,l)^ Vaer'^. (2.11) 

J a 

Furthermore, it follows from (12. 8p . (12. 7p and (I2.10p for any r > 1 and any 
a eT'^ that 

I / \tvr - (ii'i^r, i)'i <CK w\ i[/^]ni,oo,. < ctk \a\ mi^^^ 

\Jve [W^'^{a)Y. (2.12) 

In addition, let = to < ^i < • • • < ^n-i < t^ = T he a partitioning of 
[0,T] into possibly variable time steps r„ := tn — tn-i, n = 1 ^ N. We set 
T := max„=i_^jv Tn and introduce 

r(-) := 1 /"/(., t) dt G L2(fi) n=l^N. (2.13) 



n Jt„-i 



We note that 



N T 

E^"l/"lo,s,n< / \f\o,s,ndt foranysG[l, 

n=l -^0 



Finally, on setting 



w'/ = P^'ln^'w^ 



OJ 5 



2] . (2.14) 



(2.15) 



we introduce M^ : 5'' — )• S*^ approximating M^ : C{Q) — )■ C{Q), defined by 
(II. lip , for any a G T^ as 



M^ivn ■■= < 



( ko 

kii + (ko - ki::, 



h s.h 



l]"" G [u7o ' ^0^ + ^] 



£,/l 



^ fci „ := max(A;o, iVvr'^wg |<^ |) 






(2.16) 
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We note that M^ is also well-defined on S^ with M^ : S'^ -^ L°°(f2), and we 
have the following result. 

Lemma 2.1 For any r]'^ G S^, we have that 

(2.17) 

Proof. It is convenient to rewrite fll.lip and fl2.16p for any rj'^ G S^ and 
for a.e. x G fi as 

Me{v'")ix) = ko+ f^li^hlh^ niin( maxK(x) + e - r]''{x),0), e) , 



(2.18a) 

M^{7]''){x) = ko+( ^'' ^-^ ~ ^° I min( max(w;^''^(x) + e - r7'^(a;), 0), e) ; 

(2.18b) 
where 

^0 < M^{r]''){x) < kl'''{x) := max(fco, |V7r'^H(x)|) for a.e. x G fi . (2.19) 

Since 

I niin(niax(a, 0),e) — niin(niax(6, 0),e)| < |a — 6| 
and |min(max(a,0),£)| < £ V a, 6 G M, (2.20) 

it follows from ( l218al b). (I2A9|) . ( 12J[5|) . (I2:4il) and Assumption (Al) that 
\M,ir]'^)-M^i7^'^)\o,oo,n 

, '^1,00 '^0 I J s.hi I I IV7 ei irv /i ei I 

< Fo - ^0^ lo,oo,n + I lY^ol - lYvr Wq\ |o,oo,Q 

< Cie-') [ |(/ - P")w7g|o,oo,f^ + ||(/ - 7r^X||i,oo,c] ; (2.21) 
and hence the desired result (I2.17p . D 
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2.1 Approximation (Q^' ) 



Our first fully practical finite element approximation is: 

m-l N h 



{Ct/) For n = 1 ^ AT, find W^ E U^} and Q" G S^ such that 



(2.22a) 
(2.22b) 



where W^ = tt^Wq. 



For any x^ ^ f^o^? we introduce the closed convex non-empty set 

K^{X^) ■= {r/^ G t/o'' : \^v'"\ < M^{P^x^) a.e. in Vt] . (2.23) 



In Theorem 12.21 below, we will show that (Q^'^), (I2.22al b). is equivalent 
to (P^'^) and (M^'^). The former is the approximation of the primal quasi- 
variational inequality: 
(P^'") Forn = 1 ^ AT, find Wl G K^{Wl) such that 

fw"^ — iy"~^ \ ^ 

{^ ^,7]^-Wl\ >{r,r]^-Wl) ^ri^eK^iWl), (2.24) 

where IV4 = tt'^Wq- 

The latter, having obtained {Wj^}^^]^ from (P^'^), is the minimization 
problem: 
(M^'") Forn = 1 ^ A^, find Q^ G Z^''^ such that 

{M^{p^wX), |q;;|) < {m^{p^wX), \v!^\) v / g z^'" , (2.25) 

where 

r /'W^ — W'"'"^ \ ^ 

Vj^'^G t/oH- (2-26) 

) oi S!", it follows that th 
n = 1 ^ N, is non-empty. 



As Vt^o is a strict subset of S_^, it follows that the afiine manifold Z'^'", 
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We consider the following regularization of (Q^'^) for a given r > 1: 
(Q^-;) FoTn = l^N, find W^, e Ul} and Q^^^eS^ such that 

"Tn J ' 

(2.27a) 



(2.27b) 



{M':{P'^W2,r) \Ql/-'Ql,,vn + (ViyZ.,/) = V / G 5' 



where W^ ^ = ti'^Wq. 

Associated with (Q^'^) is the corresponding approximation of a gener- 
alised p-Laplacian problem for p > 1, where, here and throughout the paper, 



r p 

(P^';) Forn = 1 ^ A^, find W^^r ^ K such that 

= {r,v'') ^V'^U^, (2.28) 

where W^ ^ = ti^Wq. 

Theorem 2.1 Let the Assumptions (Al) and (A2) hold. Then for all r e 
(1,2), for all regular partitionings T^ of Q, and for all Tn > 0, there exists 
a solution, W^^r ^ ^o" ^^^ Q"]^^^ S^ to the n^^ step of ((t^,^)- In addition, 
we have that 



N N 



n=l n=l 



max^ \W2M + Y. \^lr - ^ZMn + E ^« I^IJS,.,f. < ^ , (2.29a) 

n=l 
N \3 

Y.^n\VWlMA <C (2.29b) 



\n=l 



where I + ^ = 1. Moreover, (Q'^X'^J, ^2.27a\ b), is equivalent to (PaI), ^KW- 
Proof. It follows immediately from fl2.27bp that 
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(2.30) 



Substituting this expression for Q^,, into fl2.27ap yields fl2.28p . Hence (P^'p), 

with f l2.30p . is equivalent to (Q^'^)- 

We now apply the Brouwer fixed point theorem to prove existence of a 
solution to (P^'p), and therefore to (Q^'^)- Let F^ : Uq -^ Uq be such that 

for any ip'^ G Uq, F^ip^ G Uq solves 

' ph h _ y^n-l \ ^ 



5 

Tn 



= (r,^') Vr/'^Gf/o'- (2.31) 

The well-posedness of the mapping F^ follows from noting that f l2.3ip is 
the Euler-Lagrange system associated with the strictly convex minimization 
problem: 

minEp^'"(r/'^), (2.32a) 

where E^'" : f/j^* — )■ M is defined by 

^,'''^(^') := ^ 1^' - WZ'W + - /[Mi^(PV'^)]-(P-^) lYr^'^rdx - (r,r/'^) ; 

(2.32b) 

that is, there exists a unique element (F^ip^) G f/^ solving (I2.3ip . It follows 
immediately from (I2.32al b) that 

(2.33) 

It is easily deduced from (I2.33P and (I2.9ap that 

FV' G ^7 ■■= iv' e f/o' : \v%,n < 7} , (2.34) 

where 7 G M>o depends on |M^^7^lo,n, |/"|o,n and r„. Hence F'^ : B^ ^- B^. 
In addition, it is easily verified that the mapping F^ is continuous, as M^ : 
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S'^ — i- S'^ is continuous. Therefore, the Brouwer fixed point theorem yields 
that the mapping F^ has at least one fixed point in S^. Hence, there exists 
a solution to (P^'^p), (ESHD, and therefore to (Q^fJ, ( K2TE\ h). 
It follows from (Km and ( 12191) that forn = 1 ^ A^ 

where, on noting fl2.19p . (]2.2bp and Assumption (Al), 



,00 



e.h/ 



max kl' (x) < C . 



(2.35) 



(2.36) 



Choosing rj = W^^ v = Q", in (I2.27al b). combining and noting the simple 
identity 

(a-b)a= -[a^ + {a- hf - 6^] M a,he^, (2.37) 



we obtain for ra = 1 — )■ A^, on applying a Young's inequality and (ll.lSp . that 
for all (5 > 

\WlM + \Wl^ - W2-'\l + 2r„ (M^iP^Wl^), \Ql/) 

Tn-\\2 



1 

r 



<\w2-'\1 + 2t, 



rn \p 



Z^''\no,r,n + -^'\^A,rio,p,n 



<\WZ'\1 + 2t„ 



p 
1 



rn IP 



S~' \r\o,r,n + -[SC.m'\YW2,rio,p,n 



(2.38) 



It follows on summing (I2.38P from n = 1 to m, with 6 = l/(C^(fi) [A;i'oo]'')) 
and noting (12351) and (121361) that for m = 1 ^ A^ 



n=l 

m. 



ra=l 



n=l 



The desired results ^M^h) follow immediately from ( 1239|) . dZSa]), (l27[i|) . 
(EinD, (1^:^ and (ESSD- □ 
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Theorem 2.2 Let the Assumptions (Al) and (A2) hold. Then for all regu- 
lar partitioning s T^ ofVt, and for allTn > 0, there exists a solution, W^ € Uq 
and Q" G S_ to the n^^ step of (Q^ ). In addition, we have that 



N N 

|iy"ln,o, + ' 

n=0->-Af 



max \Wl\o,n + Y.\^A- Wl-'\ln + E ^« \Ql\^X^ ^ ^ ' (2.40a) 



71=1 n=l 



max ||iy^||i,oo,Q<C. (2.40b) 



Moreover, (Q)C), ^2.22a\. h). is equivalent to (F^/), [2J4\ ), and (M^/)^ 



\2.25\) . Furthermore, for n = 1 ^ N, having obtained W^, then Q^ = 
— A^ VW^^, where A^ G S^q is the Lagrange multiplier associated with the 
gradient inequality constraint in (P^ ). 

Proof. It follows immediately from fl2.29ap . on noting that | ■ |o,i,q < 
|f2| + I ■ Ioj. ^5 that for fixed T^ and {Tn}n=ii there exists for n = 1 — ;■ A^ a 
subsequence of {l^I,„Q^ Jr>i (not indicated) and W^ e U^ and Q^^ G S^^ 
such that 

Wl^^Wl, Q^r^Ql asr->l. (2.41) 

We now need to estabhsh that {W^^,Q^}^=i solves (Q^'^), fl2:22il bl Noting 
02.4ip . one can pass to the limit r — )■ 1 in (I2.27ap to obtain 02.22ap . Choosing 



^" = 91r~^ '^^ ( IMTbjl and noting fITTU]) . fITTg]) and ([23SD, one obtains 
that 

> iM^iP'w2,r), \QIJ - 1 \t'n + ^ C 1^1 

V ^^ G 5^ . (2.42) 



Noting (EH]), ( I2.18bp and f lOOj) . one can pass to the limit r -^ 1 in (jMl]) 
to obtain fl2.22bp . Hence there exists a solution to (Q^'^), (I2.22a[ b). 



In addition, one can pass to limit r — )■ 1 in (I2.29ap . on noting (I2.4ip . to 
obtain the desired result (I2.40ap . 
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Choosing v'' = and 2Q^ in ( 12.22bl) . yields ioi n = 1 ^ N that 

(^h(ph^n)^ |gn|) ^ (Vl^^^Q^) = Q (2.43a) 

and hence that {M^{P''W2),\v!'\) + iYW^^v'') >0 V / G 5\ 

(2.43b) 

Choosing 



v^ = { —'^ '"* - -*, ^2.44) 



'-k 1 



in fl2.43bp . and repeating for all o"^ G T^, yields for n = 1 ^ N that 

IVW^I < M^iP'^W^) a.e. on f] . (2.45) 



As W2 G f/o^, it follows from ([235]), fl2:T9|) . (I2:36|) . (iLTSj) . (l2:2b|) and our 
choice of VTjJ that the desired result f l2.40bp holds. 

It follows from fl2:i5D that W2 G K^{W2). Choosing r]^ = ^^^ - Wl for 
any ^^^ G is:^(W2) in (ESSE]), we obtain, on noting ^M^ and ([223]), that 



h 



= (q:;,y(/ - wx)) = {M^{p'^w2), |g;;i) + (y/,q;;) > o. (2.46) 

Hence {W^j^^i solves (P^'"), ([221. It follows from ([2221]) that Q^^ G Z'^'", 
n = 1 — 7- A^. Therefore fl2.22bp immediately yields (I2.25|) on choosing 
^h ^ g^,n^ gg^^g {Q''Jn=i solves (M^), (1^21]) . Therefore a solution 

{W^I,Q:;}1i of (Q:^'-) solves {F'/) and (M;;,'-). 

We now prove the reverse. If {W2}n=i solves (P^'""), then, forn = 1 — )■ A^, 
W^ is the unique solution to the strictly convex minimization problem: 

min E^^'^iri^) , (2.47a) 

where E^^'^ : [/q — )■ M is defined by 

^'•"(^') := ^ 1^' - Wl-\ - (r, r/'^) . (2.47b) 
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Next we introduce the Lagrangian L'*'" : Uq x S'>q — > M defined by 

lh,n^^h^ /) := E'^'^'ir]'') + i (/, I Vr^'^p - [M,'^(P'^iy^)]2) . (2.48) 

As fco > 0, we note that the Slater constraint quahfication hypothesis, see 
e.g. (5.34) on page 69 in Ekeland and Temam [9], is obviously satisfied; that 
is, there exists an rj^ G Uq such that IVt^qI < M^{P'^W2) ■ Hence it follows 
from the Kuhn-Tucker theorem, see e.g. Theorem 5.2 on page 69 in Ekeland 
and Temam [9], that there exists a A^ G S'>q such that 

L^'^'iW^, /) < L^'"(Vr^, A;^) < L'^'"(r7\ A;^) ^7]^ e Ul^, V / G 5^0 • 

(2.49) 

The first inequality in (I2.49P yields for /x'^ = and 2A^ that 

(Al, \VW2\^ - [M^iP''W2)f) = 
=> {XI \VW2\, \VW2\ - M^{P''W2)) = , (2.50) 

as W2 e K'^iWi). The second inequality in (ITi^ yields that 

(^— ^— -^,r/M + {X\VWI,Yj1^) = {f\^'') yr^'eU^. (2.51) 

It follows that <K22^ holds on setting Q^ = -X^VW^, and Q"^ G Z'^'". 
Furthermore, we have from this definition for Q". G Z^'"' and (I2.50p that for 
all v!" G Z'^'" ~ 

(Mi^(P^W^^), \Q^\) = -{Ql,VW2) = -{v\YW2) < {M^{P'W2), \v'\) , 

(2.52) 

where we have recalled that W2 G K^{W2) for the last inequality. Hence, 
for n = 1 — )• A^, Q" = —X'%V_Wa ^ ^^'" solves the minimization problem 
(M^'"), ([22SD. Since the inequality in (E32D hol ds for a ll v!" G S'^, it follows 
from this and the first equality in (I2.52p that (]2.22bp holds. Therefore a 
solution {iy^,Q"}^^i of (P^'") and (M^'") solves (Q^'"). D 



ih,T^ 



2.2 Approximation (Q^' ) 

Our second fully practical finite element approximation is: 
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(Q^'") Forn = l^N, find WJ^ G S'' and Q^ G Z'' such that 

Tn J 

(2.53a) 

(Mi^(iy^), |t;'^| - ig^l)^ - {W^,V. (v'^-Ql)) > V/ G Z^ 

(2.53b) 

where W^^ = P'*[7r'*w^]. 

For computational and theoretical purposes, it is convenient to consider 
the following regularization of (Q^'^) for a given r > 1: 
(Q^-;) Forn = 1 ^ A^, find W^^^ G S^ and Ql^eV^ such that 

- PFS"^ \ 

(2.54a) 

(2.54b) 

where W^^^ = P'^[7r'^wg]. 

Theorem 2.3 Let the Assumptions (Al) and (A2) hold. Then for all r G 
(1,2), for all regular partitionings T^ of Q, and for all Tn > 0, there exists 
a solution, W^^. G S'^' and Q" G V^ to the n^^ step of (Qb^j.)- In addition, 

we have for r G (0, |] that 



N N 

liySJno + ' 
n=0->-N 



max^ \WlXn + Yl \^B,r ' W^-'Hn + E ^" \Ql/o,r,n < C , (2.55a) 

n=l 

N 

T.^n\Y-Qlr\ln<C- (2.55b) 



n=l n=l 

N 



n=l 



Proof It follows from (I2.54ap and 023]) that 

y^B,r = 9^-rnY-Ql,, where g^ = WJ^-' + t, P"" T . (2.56) 
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Substituting fl2.56p into f l2.54bp yields that the n*^ step of (Q^'^) can be 
rewritten as find Q" G V^ such that 



■B,r 
hf„n _ v7 /O" N l/^n \r—2/^n „.h\h i „ /'vj ^n vj „,h\ 



(M,^(^" - r„ V . g^ J IQ^ J'-^g^^^, i; Y + r. (V . Q^^^, V . i; 



= (^",Y.t;^) Vt;'^GZ. (2.57) 

One can apply the Brouwer fixed point theorem to prove existence of a so- 
lution to (12371) . and therefore to (Q^'^. Let G^ -.V!" ^ V^ be such that for 
any ^^ eV!", G^^ e V!" solves 



(M,'^(^" - r„ V .<) |G>^r-"GV,r) + r„ (V. (GV), V.t;'^) 

= (^",Y.i;^) Vi;^eV:''. (2.58) 



The well-posedness of the mapping G^ follows from noting that fl2.58p is 
the Euler-Lagrange system associated with the strictly convex minimization 
problem: 

min Jl!'''{v!') , (2.59a) 

where J^'" :]/''—)■ R is defined by 

JM.(^^) := 1 (M^(gn - r„ V .^^), |/r)' + ^ \Y.v%n ' (^", Y./); 
r — 2 

(2.59b) 

that is, there exists a unique element {G^ip^) G V'^ solving fl2.58p . It follows 
immediately from ^M^h) that J^'^'"(GV'') < J^"(0), and this yields, on 
noting fl21&|l and fl2llD . that 

^ ( |GVr, 1) + y lY. (GV)l^,f. < (^", Y. (GV) ) • (2.60) 

It is easily from fl2.60p that 

GV e 5^ := {/ G Z' : |^!'|o,r,n < 7} , (2.61) 

where 7 G M>o depends on |5'"|o,c, ^ and r„. Hence G^ : 5^ — > 5^. In 
addition, it is easily verified that the mapping G^ is continuous, as M^ : 
S^ — 7- S^ is continuous. Therefore, the Brouwer fixed point theorem yields 
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that the mapping G^ has at least one fixed point in B_^. Hence, there exists 
a solution to (Q^'^), fl2.54a[ b). 

Choosing t]'^ = WJ^^, v'^ = Q" in fl2.54a[ b). combining and noting fl2.37p 
yields, similarly to (12.38 p . that 



\W] 



n |2 
B,r\0,n 



+ Wh. - W-^-}\in + 2r„ {M^{Wl,), \QIJ) 



■\n \r\h 



\W'^:X,n + 2r„ (r, Wl,) < \W-;X^ + r„ [\WlM^ + \fX,n] 



(2.62) 



It follows from flTB2]) . on noting that (1 - x„)-i < (1 + 2r„) < e 
r„ G (0, i] and (EH]), that forn = 1 ^ A^ 

IW" P < p2r„ r |Ti/n-l|2 _, |i-n|2 1 



as 



< e^*" 



N 
|TT/0 |2 I "^T- |fm|2 

r?i=l 



<c. 



(2.63) 



which yields the first bound in (I2.55ap . Summing (I2.62p from ra = 1 — )■ iV 
yields, on noting (I2.19p . (12. lip and (I2.63p . the second and third bounds in 
( I235al) . 

Choosing 77'^ = V . Q'% ^ in (I2.54ap yields that 



'-B,r' 






Summing (I2.64p from n = 1 — )■ A^, and noting (I2.55ap and (I2.14p . yields the 
desired result ( I2.55bp . D 

Theorem 2.4 Let the Assumptions (Al) and (A2) hold. Then for all regu- 
lar partitionings T^ ofQ, and for alWn > 0, there exists a solution, W^ G S^ 
and Q" G V^ to the n^^ step of (CfB^ )■ In addition, we have for r G (0, |] 
that 



N 



N 



n=0->-N 



71=1 



maxjW^\,,n + J2\^b- W^Xn + E^" \9.>xn <C. (2.65a) 

n=l 

N 

Y.<\'^-Tb\I^^C- (2-65b) 



n=l 



22 



Proof. Similarly to fl2.4ip . on noting that | ■ |o,i,n < 1^1 + I ' lo rQ; i^ follows 
from fl2.55a[ b). that for fixed T^ and {Tn}n=i, there exists a subsequence of 
{W^B,r, Ql.r}r>i (not indicated) and W^ E S'' and Q^ e V^ such that 

WS,r^WS, Ql^->Ql asr^l, (2.66) 

and the bounds fl2.65a[ b) hold. One can now immediately pass to the limit 
r -)■ 1 in f l2.54ap to obtain f l2.53ap . Similarly to (12.42 p . choosing v^ = 



Ql - i)'^ in fl234b|) and noting f l230|) . fl2:T9|) and flOej) . one obtains that 



V ^^ e Z'' • (2.67) 

Noting (12. 66 p . one can pass to the limit r — )■ 1 in (I2.67p to obtain (I2.53bp . 
Hence there exists a solution to (Q^'^), (I2.53al b). D 



3 Convergence 

We introduce the following discrete approximation of the mixed formulation: 
(Q^) Forn = 1 ^ AT, find w" G WI'°^{VL) and f G V_^{Vl) such that 



Tn 



r? +(V.g",r/) = (r,r7) Vr7eL^(fi), 



(3.1a) 

(Il!| - |g1,M,(^"))c(^) - (V. {v-f),w-) > V t; G V^iSl) ■ 

(3.1b) 

where m;° = iUq. 

For any x ^ W^o'°°(^)' ^^ introduce the closed convex non-empty set 

i^(x) := {r/ G ^■'"(fi) : iVr^l < M,(x) a.e. on fi} . (3.2) 

Associated with (Q"^) is the corresponding approximation of the primal quasi- 

variational inequality: 

(P^) Forn = 1 -^ A^, find w" G K{w'^) such that 

^"~^" ,r]-wA > (r,r/-w") V r/ G ir(ti;") , (3.3) 
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where w'^ = w^ 



0- 



In Section [STT] we show, for a fixed time partition {t„}^^^, that a subse- 
quence of {{W2,q'X=,}h>o, where {W2,Ql}^=, solves (Q^^^'^), converges, 
as /i — 7- to {w",g"}^^i solving (Q"^). In Section [3^ we show, for a fixed 
time partition {r„}^^^, that a subsequence of {{W^s^,(5" }n=i}h>o, where 



{^B,r,QlJn=i solves (Q^-;), couvcrgcs, as /i ^ and r ^ 1, to {w'^, g"}^^i 
solving (Q^). For our final convergence result in Section 1231 we need an extra 
assumption on the data. 

(A3) u;g > and / e L^{0,T; L^n)). 

Under this further assumption, we will show that a subsequence of 
{{w^,q"'}^^i}r>o, where {w",5'"}^=i solves (Q"^), converges, as r — )■ 0, to 
{w, q} solving 

(Q) Find w e L°°(0,T;iyo'°°(^)) ^ iy^'°°(0,T; [C^m]*) and q G L°°(0,T; 
[Min)Y) such that 



,dw 



v) 



dt ' '/^o(^) 



hY.v) 



cin) 



if,v) 



dt = Q 

\/r]eL\0,T;C',m), (3.4a) 



{\v\-\q\,Me{w))c^^~{V.v-f,w) 



dt 



>UK^T)\l,n-Hi-)\ln] yveL\0,T;V^m; (3.4b) 

where w{-, 0) = Wq{-). 

Associated with (Q) is the corresponding approximation of the primal 
quasi- variational inequality: 
(P) Find w G L~(0,T; Kiw)) H W^'^{0,T; [C^m]*) such that 

CJIV 

\-g^'V)c^^(n)-if,V-w) dt> I [k(-'^)lo,n-ko(-)lo,f7] 

\/ r]eL\0,T;K{w)r\C^m), (3.5) 
where w{-, 0) = u^o(')- 

Remark 3.1 One might expect the inequality in the primal quasi-variational 
inequality (P) to he such that 



T 



dw 



{f,v-w) 



dt>0. 



(3.6) 
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However, the term 

''^ ,dw 



^.^)ci(fl)dt (3.7) 

IS not well-defined for w E L°°(0,T; Wo'°°(f]))nM/^'°°(0,T; [C^m]*), and has 
been rewritten to yield /i3. 5\) . which is well defined. This follows from (E 



with B = L'^{fi), and, for example, the reflexive Banach spaces Bq = Hq{Q) 
and Bi = [VFo'^(^)]* '^'^^^ s G (d, oo); see the proof of Theorem \3.5\ below. 
In addition, the test space has been smoothed to make the first term on the 
left-hand side of Ii3.5\) well-defined. 

Similar remarks apply to (Q), where one might expect the inequality in 
( 3.4b ) to take the form 



I \{\v\-\q\,M,{w))c(^)-{V.{v-q),w)\dt>^. (3.8) 

Jo '- -■ 



However, the term 



[ iY-q,w)dt = - [ (g,Vw7)c(n)dt (3.9) 

Jo Jo 

IS not well-defined for w G L°°(0,T; WQ'°°{n)) and q G L°°(0,T; [Mm^]. 
This term has been rewritten using Iji3.4a\ ) formally with rj = w, and the 
rewrite of the term ( [5'. 7| ) employed in Ii3.5\) . to yield Ii3.4b^ . 



3.1 Convergence of (Q^' ) to (Q^) 

Theorem 3.1 Let the Assumptions (Al) and (A2) hold. For any fixed time 
partition {t„}^^;^ and for all regular partitionings T^ of Vt, there exists a 
subsequence o/KVF^, (5'!}^=i}h>o (not indicated) , where {W2,Q'^}^^i solves 
(Qa)' '5'wc/i that as h —^ 

VW2 -^ Vw" weak* m [i:~(fi)]^ 

W2 -^ w" strongly in C(fi), 

M^{P''W2) -^ M.iw'') strongly m L~(fi), 

qa _^ ^n vagucly in [M{Vt)] 



d 



n = 0^ N, 


(3.10a) 


n = -> iV, 


(3.10b) 


n = 0^ N, 


(3.10c) 


n = l^ N; 


(3.10d) 



where {w"", 5'"}^=i is a solution of (Q'), Ii3.1a\ b). 
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Proof. The desired subsequence convergence results fl3.10al b.d) for a fixed 
time partition {t„}^^^ follow immediately from fl2.40a[ b). on noting that 
W^'°°{Q) is compactly embedded in C{Q) and fll.lSp . Next we note that 

< iM^iP'^W^) - M,{P^W2)\o,oo,n + \M,{P''W2) - M,(u.")|o,oo,q • 

(3.11) 

It follows from ( I2.18ap . (^M) and fl^:ial) that 

<C{e-')\w^-P''W2\o,oo,u 

< C{e-') [\{I- P")«;"|o,oo,Q + k" - W2\o,oo,n] • (3.12) 

Hence, the desired result (IXTU^ follows from (IXTTj) . (IZT7|) . (IXT^ . (E^b]), 
f l2:2b|) and flSlObl) . 

We now need to establish that {if", g*^}^^]^ solve (Q'^), (I3.1a[ b). For any 
rj G C^{Q), we choose r^'^ = tt^t] in (I2.22al) and now pass to the limit /i — > 
for the subsequence to obtain, on noting (13.10bp . (]2.9bp . fl2.2bp and (IS.lOdp . 
for n = 1 — )■ A^ that 

'^\v)- {f,Y.v)c(n) = ir,v) "^Ve C^m. (3.13) 



w — w 



Tr,. 



It follows from ( IXTSj) . (l2A4l) and as w" G C{VL) that 



(f,Yr/)^(^) <C(r-^)l^ 



o,n 



^ Tie C^{Q) 



(3.14) 



We deduce from 03.141) that the distributional divergence of g" belongs to 
L^(r2), and hence g" G I£"^(fi), n = 1 — )■ iV, and so (I3.13p can be rewritten 
as 

,n-l 



W — W 



r/ +(v•g^^) = (/^^) vr^GCo^(fi). 



(3.15) 



Noting that C^{n) is dense in L'^{Q) and that tt;", V .g'^, /" G ^^(fi) yields 
the desired 03.1al) . 

For any v G [C°°(r2)]'^, we choose v^' = P_^v in (12. 2 2b P and now try to 
pass to the limit for the subsequence as /i — )■ 0. First we note from (I3.10ap . 
f l2:4b|) and as w" G 1^0'°° (fi) that forn = 1 ^ iV 

\im{VW2,P^v) = (Y«;",t;) = -(w", V.t;) . (3.16) 

/i->-0 
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It follows from (l^^^i]) with r]'' = W^, (l3J0b|) . (l2Jb]), I^M^ and ( KTa^ 
with r] = w" that for n = 1 ^ N 



lim(VVr^,Q") = lim 



h^O 



Next we note that 



h~>0 






w — w 



Tn 



n-1 



Tn 



,W 



{FM^ 



[w\V.q^). (3.17) 



{M^{p^wi), |q;;i - \p^v\ 



(M,(^"), |Q"| - \PW) + (M,'^(P'^iy^) - M,(u;"), |Q"| - |PS|) . 



lilyll 



(3.18) 



As M£(tf;") G C(fi) is positive, it follows from (13.indl) . (ILT9D and ( l2lbll that 
hminf (M,(t."), |Q;;| - \P\\) > {\f\ - \v\,M,{wn)cm ■ (3.19) 

It follows from (l2.4Ua|) and (EH]) that 

(M^'^(P'^Vr^) - M,(w;"), |Q;;| - iP'^ul) 

< \M';{P''W2) - M,(i/;")|o,oo,f7 [Cr-i + Ho,i,n] • (3.20) 

Combining (I3.16p - p.20p and (I3.10cp . we can pass to the limit for the subse- 
quence as /i -;■ in (I2.22b|) . with v^ = P^v for any fixed v G [C°°(Ji)Y, to 
obtain for n = 1 ^ N that 

{^.{f-v).w-) > (If I - \vlM,{wn)cin) ^velC^iQ)]". (3.21) 

Recalling the density results (ll.21al -c) and that w", M^{w"') G C{il), we 
obtain the desired result ( I3.1bl) . D 

3.2 Convergence of (Q^^ to (Q^) 

For the purposes of the convergence analysis in this subsection, it is conve- 
nient to introduce the following regularization of (Q"^) for a given r > 1: 
(Q;) Forn = 1 ^ A^, find < G W^'^in) and g" G V'in) such that 



w" — w" ^ 
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(3.22a) 



(M^iw:) \q:r\:,v) - «, v .1;) = yve rm ■, (3.22b) 



where w^ = Wq. 



Theorem 3.2 Let the Assumptions (Al) and (A2) hold. For any fixed r e 
(1,2) and fixed time partition {Tn}n=i with r G (0,^], and for all regular 
partitionings T^ of Q, there exists a subsequence of {{WJ^j.,Q'^ }^=i}h>o 

(not indicated), where {W^^., Q^ }n=i solves (CfB^r)' such that as h -^ 0, for 

any s G [1, oo), 



^B,r ^ '^r strongly in L^(fi), 

M^iWB,r) -^ M,{wl^) strongly m L'(fi), 
Ql^-^q: weakly zn[U{n)] 



d 



-LB,r 



V.g" ^V.g" weakly mL\9), 



= 0->iV, 


(3.23a) 


= 0^ A^, 


(3.23b) 


= 1 ^ A^, 


(3.23c) 


= l^N- 


(3.23d) 



where {w'^^q^}^^^ is a solution of (Q^), \3.22a\, h). 

Proof. The desired subsequence weak convergence resuhs fl3.23cl d) follow 
immediately from the bounds on {Q^^}^^-^ in fl2.55a[ b). on noting that the 
time partition {Tn}n=i is fixed. In addition, we obtain from the first bound 
in ( I2.55al) that 

W^B,^< weakly in L2(fi), ^ = ^ A^. (3.24) 



Furthermore, we obtain from f l234b|) . ^U^ . 0231]), (EIID and f E3^ for 
n = l^ N that 

mi,.^.v^)\ = \{M^{wi,) \Qi/-'Qiyvy'\ < ki:i{\Qi/-\ ia)" 

< c mi/, 1)1^ [(l/r, 1)1^ < c ig^js;;^ |t;^|o,.,n 

< Cir/) \v%,r,n yv^eV''. (3.25) 

For any fixed v G [C°°(i7)]'^, on choosing t/" = I^v_ in ( I3.25p . letting /i — > 
and noting ([2l]), (13:241) and ([22D, we obtain that 

|«, V.t;)| < C(r-i) Ho,.,n, n = l^N. (3.26) 
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Repeating (IS:^ for all v G [C~(n)]'^ and as C°°(fi) is dense in U{Q), we 
obtain that 

w'^eW^'^'in) with |K<||i,p,Q<C(r„"i), ^^i^jv, (3,27) 

The fact that w" vanishes on 5^ can be deduced from fl3.26p by using an 
argument similar to that in [6l p. 699]. 

Next, for ra = 1 -)■ iV, we introduce V^W^^,r ^ I^^ such that 

(Y^^^,„i;^) = -(Vr^,„V.i;'^) ^v'^eV''. (3.28) 

It follows from ( K2E\> and (ESSD that 

|Y^W^^,.|o,Q < C(r-i), n = 1 -> AT . (3.29) 

For n = 1 -^ N, we now introduce l^g r ^ ^o such that 

(Y%,, Vr/'^) = (V,iy^,,, Vr/^) V r^^ e f/o' . (3.30) 

It follows from ( ITTB]) . (1^:^?]) and dSSHD that forn = 1 ^ A^ 

||W^B,rlli,n < C7|W5^,|o,f7 <C\Vf,Wl,\o,n < C(r-i) . (3.31) 

We deduce from fl3.29p and f l3.3ip that there exists a further subsequence 
of {{V_h^B ry^B r}n=i}h>o {'^ot indicated) such that as /i — )■ 0, for any 

s e [i,oo), 



YhWI^ ^ € 


weakly in [^^(fi)]'^, 


ra = 1 ^ AT, 


(3.32a) 


Yw;^^, -^ YWr 


weakly in [^^(fi)]'^, 


n = 1 ^ A^, 


(3.32b) 


Wl, ^ W^r 


strongly in L*(fi), 


n = 1 ^ AT; 


(3.32c) 



where w" G //({(i^). For any fixed y_ G [C°^(r2)]'^, on choosing t;'^ = / t; in 
(13:281) . letting /i ^ for the subsequence and noting (ESI), ( l3:32all . (JSSID 
and (12. 7p yields that 

(^",^) = -«,Y.w) n = l^Ar. (3.33) 

Repeating ( 1333|) for all t; G [C°°(n)]'^ yields that d;? = Vw;?. Similarly, for 
any fixed ?7 G C^(r2), on choosing r^'^ = i^^ri in f l3.30p . letting /i — )■ for the 
subsequence and noting (]2.2bp . f l3.32a[ b) and ^ = Viu" yields that 



{Ywr, Y.V) = (^", Yv) = (Vw", Vj]) n = l^N. (3.34) 
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Repeating ( IXM]) for all rj G C^(fi) yields that w^^ = w"^. 
For n = 1 — )• A^, let z" be such that 

-Az" = W^^, - Wl^^ in f], z" = on 9f]. (3.35) 

As n is convex polygonal, elliptic regularity yields that 

lk"||2,c<C|l?^,,-iy^,,|o,c. (3.36) 

It follows from ([S35D, dSSOD, (ESD, (IS2HD, dXM]) . (1^^ . (1^ and ^M 
that for n = 1 — )■ A^ 

= (Yvr^,.,Y^") + (iy^,„A^") 

= (W^,., Y(^" - ^z-) ) + (V^iy^,,, V[7r'^^^) + (ly^,,, Az") 

= ( W^,. - Y,iy^,„ V(z" - vr^^") ) + (V,W^^,„ V^") + (W^^,,, A^") 

= (Yl?^,, - Y.W^5,., Y(^" - vr'^") ) + (Y.W^^,., Y^" - I'(Y^") ) 

< C(r-i) [ |V(z" - 7r^^")|o,n + lY^" - l'(Y^")|o,n] 

< C(r-i) h \\z\\2,^ < Cir^') h\ (3.37) 

As w" = w", ra = 1 — )■ A^, it follows from (I3.37P and (I3.32cp that the desired 
result (13:2^ holds. 

We deduce from (I3.23ap . (I2.18ap and (I2.20p for a further subsequence of 
{{^B,r]n=o}h>Q (not indicated) that as /i -)■ 0, for n = -)■ A^, 

ly^^, ^< a.e. infi ^ M,(Vr]^,,) ^ M,«) a.e. in fi . (3.38) 

It follows from (13.380 . (I1.14p . (I1.13P and the Lebesgue's general convergence 
theorem that as /i — )■ for any s G [1, oo) 

M.iW"^^) -^ M.iw"^) strongly in L'(fi), n = Q -^ N . (3.39) 



Combining (ETTp . (ITibl) . (l2:2b|) and ( IX^ yields the desired result ^^^. 



We now need to estabhsh that {w", g^j^^i solve (QJ^), (I3.22a[ b). For any 
7] G C|^(f2), we choose r]^ = P^f] in (I2.54ap and now pass to the limit /i — > 
for the subsequence, on noting (I3.23a[ d) and (]2.4bp . to obtain (I3.22ap for all 
ri G C^in). Noting that C^i^) is dense in L'^{n) and that <, V.g^^, /" G 
L2(fi) yields the desired fl5^^ . 
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For any v G [C°°{Q)Y, we choose v'^ = Q^l - I'^v in (12.54bp and now try 



■B,r 



to pass to the hmit for the subsequence as /i — )■ 0. First, we note from f l2.10p 
and fl2JT]) that ior n = 1 -)■ N 



>^iM^iWl,),\Ql/~\tvn'' 

>l:iM^iwi^)AQl/-\tiin 

+ ^ [{M^{Wl,),\tvn-{M^{Wl,.),\tvn''] ■ 

(3.40) 

Once again, it follows from f l2.10p that 

(3.41) 
In addition, it follows from ( l219ll . ^HM and ( l212ll that 

J \{M^{Wl,,)AL\Y)-{M^{Wl,)AL\Yf\<Ch\\y_\\,,^^^. (3.42) 

Combining (I3.40p and (I3.4ip . and passing to the limit /i — ;■ for the sub- 
sequence yields, on noting (12. 6p . (I3.23al -d). (12. 7p and (I3.42p . yields for n = 
1 -^ N that 



K, v. {q:-v_) > {M,{w:) \vr'v,q:-v) Wve [C^{n)f . (3.43) 



As <, Me«) e C{Q), g^ e Z''(fi) and /" G ^^(fi), it follows from (OH]) 
that (I3.43P holds true for all v E VJi^fl). For any fixed z e V_^{^), choosing 
V = q"^ ± az with a G M>o in (I3.43P and letting a — )■ yields the desired 
result ( I3.22bp on repeating the above for any z G VJi^)- Hence {if", g"}^=i 



is a solution of (Q;), dS^SUb). D 

Theorem 3.3 Lei the Assumptions (Al) and (A2) hold. For any fixed 
time partition {r„}^^^ with r G (0, |], i/iere ea;isfo a subsequence of {{w", 
g"'}^=i}r>i ('j^oi indicated), where {w^, q'"'}^=i solves (Ql), such that asr — )■ 1 

w" -)■ w" strongly in C(Jl), n = -^ N, (3.44a) 

31 



n = O^N, 


(3.44b) 


n = l^ N, 


(3.44c) 


n = l^ N; 


(3.44d) 



M^«) -^ M^(w;") strongly in C{U), 
q"" -)■ g" vaguely in [M{^)]'^, 

V-q^^Y-f weakly m L^ (fi) , 

where {w" , q"'}^^i is a solution of (Q^), Ii3.1a\ b). 

Proof. It follows immediately from (I2.55a[ b). (I2.10p and (I3.23al c.d) that 

N N 

max K|o,n + Yl \< - <"'lo,n + E ^- ICIo,.,^ < ^ ' (3-45a) 



n=0-i>iV 



n=l n.=l 

AT 



E^'l^-Clo.n<C^- (3.45b) 



n=l 



The desired convergence results (I3.44al -d) follow immediately from (I3.45a[ b) 
and (13.271) on recalling that the embedding W^'^^Q) M- (7(^2) is compact for 
p > d, Mg : C{Q) — )■ C{Q). One can immediately pass to the limit r — )■ 1 for 
the subsequence in (I3.22al) . on noting (I3.44a[ d). to obtain (13. lap . Similarly to 
dm, choosing v = q'^-tpm (^22\^ and noting fl2:T0|) . firTij) and (jSlSaD, 



— r 



one obtains for n = 1 — )■ A^ that 



«, v. (g; - V;) ) = (M,«) IgJ-^g;, g; - 3^) 

> (M,(^"), |g"|) - C(r-i) |M,(«;„) - M,«)|o,oo,f7 



- ^ (M,«), IV^D + i^ A;f_^ \n\ V V; G r (fi) • 

(3.46) 

Noting (I3.44a| -d) and (ll.lQp . one can pass to the limit r — > 1 for the subse- 
quence in IKm to obtain fl3lbl) . Hence {w'',f}^^^ solves (Q^), (]31itb). 
D 

Remark 3.2 It appears necessary to split the convergence proof of solutions 
of (Qb\) to solutions of (Q^ ), as h —> Q and r — )• 1, by first considering 
the limit h ^ to solutions of (Qr), then the limit r — )■ 1 to solutions of 
(Q^ )■ Similarly, it does not appear possible to directly prove convergence 
of solutions of (Q^^) to solutions of (Q^ ), as h ^ 0. For example, if we 
attempted the latter, we would still only be able to show Ms(WJ^) — >■ Ms{w"') 
strongly in L'^{Q) for s G [1, 00), as h ^ 0; and this is not adequate to pass 
to the limit h -^ in (Q)f), ^MMb). 
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3.3 Convergence of (Q^) to (Q) 

First we note the following result. 

Theorem 3.4 Let the Assumptions (Al) and (A2) hold. If {w"^ , q"'}^^^ is 
a solution of (Q^), ^S.ld h). then {w"}^^i solves (F'), [37^) . and 

W7">w"-^ n=l^N. (3.47) 

Proof. Similarly to f l2.43al b). we deduce on choosing v = and 2g" in 
(l3Tb|) that 

{\f\,M,{wn)cm = iY-f,wn (3.48a) 

andhencethat {\v\, M,{w''))c(n) > {'V.v,w'') Vi;GZ^(fi). (3.48b) 

Notin g that [0^(11)]'^ C V-^{n) and w"" G Wo'°°(fi) C C(n), we deduce 
from (l3^8b|) that 

(M,(ti;"), |t;|) > -(V^",i;) ^vE [C^iU)]'' , (3.49) 

It follows from f l3^ that for n = 1 -^ A^ 

|Vw"| < Me(u;") a.e. on fi ^ ^" G ir(w") ; (3.50) 

see, for example, the argument in |6, p. 698]. 

Similarly to (I2.46p . choosing t] = ip — w" for any ip G K{w'^) in (13. lap , we 
obtain, on noting (I3.48ap . employing a sequence of the form (ll.21al -c) (with 
V and {Vj}j>i replaced by g" and {q^}j>i, respectively) and (13. 2p . that 



w" — w^ ^ 



f^^-w' 



= _(V.g",y,-^") = (|g"|,M,(^"))c(^) - lim(V.g",¥.) 



J-^OO 



= (|g"|,M,(t/;"))^(^)+lim(g",V^) 



> {\q%M,{wn)c^n) " lim (|g;|, M,(ti;"))^(^) = 0. (3.51) 



Hence {w"}^=o solves (P^), (lOll . 

Let 7] = w"' + [w"'~^ — w"]+, where [s]+ := max(s, 0) for any s G R. It 
follows from (I3.50p and (I1.14p that for a.e. x G f2 

w'^ix) > w''-\x) => \Vr]{x)\ = |V«;"(x)| < M,(«;"(x)) , 
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< M,(ti;"(x)) . (3.52) 

Hence t] = w"" + [w"'~^ — w"']+ G K{w"'). Substituting this into f l3.3p . and 
recalling that the source /" > yields for n = 1 -^ N that 

I K"^ - w^U lo,n < -rn (r, K"' - ^1+) < , (3.53) 

and hence the desired result 03.47p . D 

Remark 3.3 We note that the monotonicity result ^3.4'^ for {w'^}^^^ solv- 
tng (Qn = {P-) does not hold for {W^^,,}li solving (QX^) = {F^p), 
{Wl}^^, solving (QY) = {PY), {W-^X=i ^^^^^"^ (Q^bY «^^ {W'^}n=i 
solving (QY )■ 

We introduce the following notation for t G (tn-i,^™], n = \ ^ N , 

f^'H;t) := r(-) W^{; t) := ^l^ll^ «;"(.) + ^-^LzA ^n-i(.) , 

^-.+ (.,t):=^^"(.), ^^'"(.,t):= «;""!(.), f '+(-, t) := g"(.) . (3.54) 

In addition, we write ly"^*^'^^ to mean with or without the superscripts ±. We 
note from flXM]) and f l^TT^ that 

r'+ ^ / strongly in L"(0, T; ^^(fi)) as r ^ ; (3.55) 

where s = 2 if Assumption (Al) holds, and s = oo if (A3) holds. 

Adopting the notation (I3.54p . (Q^), (I3.1a[ b). can be rewritten as: Find 
w^ G L~(0,T;Wo'~(fi))nW^i'°°(0,T;L2(f])) and g"'+ G L°°(0,T; [M(p)Y) 
such that 



/^ [(^,^) - r ^Yry)c(n) - (r^^) 



dt = 

Vr/GLi(0,T;Co^(n)), (3.56a) 



fT 



dt 



I [(N - \tYMe{w^Y)cm - {^■V-f^'^,W^' 

>|[K(-'^)lo,n-K(-)lo,f^] Vt;GLi(0,T;Z^(fi)); (3.56b) 

where w'^{-,0) = Wq{-). 
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fl3.56ap is obtained from f l3.1ap by choosing ri{-) = f^" x(-, t) dt in fl3.1ap 
and summing from n = 1 ^ N, for any x ^ -^^(0, T; Cq{Q)), and noting that 



/ (Y.f'+,x)dt = - / (f'+,Vx)c(n)dt ^x^L\0,T;C',{n)). 
Jo Jo 

(3.57) 

Similarly, (]3.56bp is obtained from (]3.1bp by choosing v{-) = ^ J^" ipi'^t) ^t 
in (I3.1bp . multiplying by r^ and summing from n = 1 -^ N, for any ip G 
L^{0,T]V^{Q)), and noting from (jSHD and (IZBTj) that 

N 

n=l 

N 

= \ [ k^lo,n - \<\ln ] + E [l 1^" - ^"~'lo,n - rn (r, w-)] 

n=l 

N 

> \ [ k^io.f^ - imo,n ] - E ^" (/^ ^'^) • (3-58) 



n=l 



Theorem 3.5 Let the Assumptions (Al), (A2) and (A3) hold. For all time 
partitions {Tn}n=i! there exists a subsequence of {{w'^,q^}^^i}r>o (not indi- 
cated), where {w", q"'}n=i solves (Q^), such that as r — )■ 

w\ w"'^ -^ w weak* m L~(0, T; W^'°°{n)), (3.59a) 

Olll^ (/111 

__ _, _ yaguely m L^{Q,T- [Clm*), (3.59b) 

w^ ^w strongly m ^([0, T]; C(n)), (3.59c) 

w^'^ -> w strongly m L'^{0, T; C{n)), (3.59d) 

M^(w^) ^ Meiw) strongly m C([0, T]; C{n)), (3.59e) 

M,(w;"'^) ^ M,(w;) strongly m ^^(0, T; C(n)), (3.59f) 

g"'+ ^ g wa^we/y m L'^{0,T; [Mm]% (3.59g) 

where {w,q} is a solution of (Q), Iji3.4a\ b). Moreover, w solves (P), /i3.5\) . 

Proof It follows from f lTTHD . f l330D . flLTJD and (Al) that 

max ||w"||i,oo,Q<C. (3.60) 
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Choosing rj = w^ in fl3.1a|] . summing from n = 1 ^- N and noting fl3.48ap . 
(1^^ . drii and IKm yields that 



N 



N 



\W 



N\2 



n=l 



n=l 



=1 



fn nN 



n=l 



Af 



2 / N 



<H\in+'2[Y.^n\niu] E 



'^n 1""^ lo,Q 



<C. (3.61) 



vn=l 



vra=l 



The bounds f l3.60p and fl3.6ip only assume the Assumptions (Al) and (A2). 
Choosing rj = w"^ in (13 .lap , and noting (I3.48ap . (I3.47p and (A3), yields 
ioT n = 1 -^ N that 



(If |,M,(t^"))^(^) = (V.f ,«;") = ^r - 
<\no,n\w^Wn<C. 
Therefore fl332|l and flTTD yield that 

max / |g"| < C . 



w^ — w'"' ^ 



T„. 



,W 



n=l^N J^ - 



(3.62) 



(3.63) 



We obtain from fl31aD . f l333D and (A3) for n = 1 ^ iV that 



w"- — w" -"^ 



r^. 



|(r-V.g",r/)|= (/^r^) + (g^Vr^) 



< 



1/ io,i,f7 + /_iq; I 



<C||^lli, 



oo,C 



'C(n) 

|^||i,oo,n 

Vr7GCo^(n). (3.64) 



Combining the bounds (I3.60p . (I3.6ip . (I3.63P and (I3.64p . we obtain, on 
adopting the notation (I3.54p . that 



I tI ±) II 



dw'^ 



dt 



L-(0,T;[Ci(n)]*) 



+ "I 11^^' — W^' IIl2(o,T;L2(Q)) + ||q;^' llL°°(0,T;[X(n)]d) 



<C. (3.65) 
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It follows from f l3.54p and the third bound in f l3.65p that 



\W 



(3.66) 



The subsequence convergence results fl3.59at b) and p.59gp follow immedi- 
ately from the bounds f l3.65p and f l3.66p . To apply f ll.23p to w"^, we first note 
that Cq{VL) is not a reflexive Banach space. However, Wq''^{VL), the closure 
of C^{VL) for the norm || ■ ||2,s,q, with s G (c?, oo) is a reflexive Banach space 
such that WQ^^iyt) C CliVt). Hence, the first two bounds in fl3.65p yield for 
s G (c?, oo) that 



\w' 



|L°°(0,T;VKl'°°(n)) 



+ 



dw'' 



dt 



<C. 



L°°{0,T;[W/"o'^(n)]*) 



(3.67) 



Next we note that the reflexive Banach space PVq '**((]) is dense in L^(i7). 
It follows that [L^(f2)]* = L'^{Q) is continuously embedded and dense in 
[VFQ''*(r2)]*; see, for example, the first two remarks in §5 in Simon J2T] . 
Furthermore, we have that C{Q) is continuously embedded and dense in 
[VFq'^I^)]*- Hence, on recalling the compact embedding of W^'^^i^Q) into C{Q) 
ioT s > d and that M, : C{n) -^ C{n), we obtain from I K67} . (IT:23|) and 
(I1.12P the strong convergence results (I3.59cl e). It follows from (ll.22p for 
s > d and a{s, d) G (0, 1) that 



\w 



w 



T,±l|2 



llL2(o,T;C(n)) 



\W — W 



r,±||2(l-a) 



lL2(0,T;L2(n))- 



(3.68) 



Therefore the strong convergence results (]3.59dt f) follow immediately from 
(ESHD, (Km> . (EUD, (1339^ and (irT2|) . 

On noting ( I3.59b[ g). Assumption (A3) and (I3.55p . we can pass to the 
limit r — 7> for the subsequences in (I3.56ap to obtain (I3.4ap . 

We now consider passing to the limit r — )■ for the subsequences in 
<K56^ . where at first we fix t; G C~(0,T; [C°^{U)]'^). Noting f l339c] -g). 
Assumption (A3) and (I3.55P we immediately obtain fl3.4bp for the fixed v G 
C°°(0,T; [C°°{Q)Y). The only term that requires some comment is the one 
involving g"^'"^, which, similarly to (I3.18p - (l3.20p . we now discuss. First we 
note that 



(|f'+|,M,(^^'+))^(^)dt 
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Jo 

(3.69) 

As Me{w) E C{[0,T],C{n)) is positive, it follows from ( [339gD and ( irTOjl 
that 

liminf / (|f'+|,M,H)^(^)dt> / (|g|,M,H)^pdt. (3.70) 

"^^•J Jo Jo 

It follows from ( Km and d339fl) that 



lim 

r-s-O 



(|f'+|,M,(u;-'+)-M,(u;))^(^)dt 

< hm C ||M,(^^'+) - M,{w)\\^.^,^T.^cm = 0- 



(3.71) 



Combining (I3.69p - (l3.7ip yields the desired result. Finally, we obtain the 
desired result (]3.4bp by noting that any v G L^{0, T; V_^{Q,)) can be approx- 
imated by a sequence {Vj}j>i, with Vj G C°°(0,T; [C°°(i7)]'^), on recalling 
(ll.21al -c): and that all the terms in (13 4bp are well-defined. Hence we have 
shown that {w,q} solves (Q), (I3.4al b). 

We now show that w solves (P), (13. 5p . Choosing w = in fl3.4bp yields 
that 

{\qlM,{w))cmdt > - 



{f,w)dt + l[\w{;T)\l^-\w^,{-)\l^] . 

(3.72) 

Then for any -q G L^(0,T; K{w) H C^iH)), on noting ( l3T2l) . we have that 



' {q,Yv)c(n)^t>- [ {\glM,{w))c(n)<it 



> 



{f\w)dt+l[\w{;T)\l^-H{-)\l^] . (3.73) 



Using the relationship (I3.73P in (I3.4ap . we obtain (13. Sp . Finally, we need to 
show that w G L°°(0,T; K{w)), as opposed to just w G L°°(0,T; iyd'°°(fi)). 
It follows from f l3^ that 



'{\q\,M,iw))cin)-if,^)] dt + | [ k(-'^)lo,n - ko(-)lo,n] 
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<J:= inf J(t;), (3.74a) 

?;eLi(0,T;y-^(C)) 



where 



J{v) 



(N,M,H)c(n)-(V.2;,w) 



dt. 



(3.74b) 



Choosing w = yields that J7 < 0. If ^7 < then, for any minimizing 
sequence {Wj}j>i, we obtain that J{2y_j) = 2J{v_j) -^ 2J < J , which 
is a contradiction. Hence JT" = 0, and so we have that J{y) > for 
any y_ G L}{Q,T]V^{Vt)). Since this is true also for —y_, and as w G 
L°°(0,T; Wl'^iVt)), we obtain that 



[v_, Vw) dt 



/ {Y-v,w)dt < / 
Jo Jo 



{\v\,M,iw))dt 



yveL\0,T;W^'\Q)); 



and therefore by a density result that 



(v, Vw) dt 



< 



{\v\,M,{w))dt 



WveL\0,T;L\n)). (3.75) 



For any p G (2, oo), choosing v = \[Me{w)]-^ Vw\p~^ [Me{w)]~^'Vw in IKTB . 
and noting the continuity of the p norm for p G [1, oo], we obtain that 

\\[Ms{w)]''^ V_W ||L-(0,T;L°°(n)) < 1- (3.76) 

Hence we have that w G L°°(0,T; K{w)), and so w solves (P), (13. 5p . D 

Remark 3.4 Under only Assumptions (Al) and (A2), we obtain in place of 
the second and fourth hounds in Ii3. 65\) that 



dw'' 



dt 



_ + ||y \\L^{0,T;[M{n)]'i) -^ ^ ■ 
Li(0,T;[Ci(n)]*) 



(3.77) 



The second bound in ^3.7'T\) follows from \3. 61\) and ^l.l^^ , whilst the first 
follows from 1^3. 64\ ) and the second bound in 1^3. 71\) . Unfortunately, the first 
bound in ^3. 64 ) is not enough to obtain strong convergence ofw'^ using (QTl 



as we require a > 1. Hence, the need for the additional Assumption (A3). 
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We believe the assumption Wq> in (A3) is not really required to prove 
l\3. 6^) . and the assumption Vwg ■K< ko on dVl in (Al) should he sufficient. 
For n = 1 ^ N, as w"' — w'^ E Co(fi) is nonnegative, it follows from (Al) 
that V_w"' .u < ko on dD.. Formally, g" = — A" Viy" in f2 with A" > 0, and 
as g" = on subcritical slopes, this yields that /^^ V . g" dx = /^^ g" . z/ ds = 
— J„^ A" Vif" .uds > 0. This can then be exploited in Ii3. 6S\) by noting that 

(|g"|,M,K))^(^) < (V.g",«;" + mt) = (f^ - '^" ~ j" \ w- + m] 

<irio,nk" + 9?t|o,f7<C, 

where VJl = max„=i^Ar Hw^^Ho.oo.q < C. Unfortunately, we are not able to 
make rigorous the formal argument above establishing that /q V.g^dx > 
under Assumption (Al). 

4 The Nonlinear Algebraic Systems 

4.1 Solution of (Q^'^) 

To solve the nonlinear algebraic system arising from (Q^'^), we recall The- 
orem 12.21 and use an extension of the splitting algorithm, ALG2, see page 
170 in Glowinski [H] from the variational to the quasivariational case. We 
introduce the Lagrangian £^'" : U^ x S_'^ x S_'^ ^- R defined by 

£M(^/^^^/^^^/^) := E^'"ir]^') - {v^Yv'" - 1^) , (4-1) 

where E^'"'{-) is defined by fl2.47bl) . For a given p G IR>o, "we introduce the 
augmented Lagrangian £^'"' : U^ x S^ x S^ ^ R defined by 

C';'-iv\t\v'') ■■= C^^''iv\t\v'') + ^ \Yr]'-t%n- (4-2) 

For any x'^ G f/g , we introduce the closed convex non-empty set 

^^(/) := {^'^ e S!" : IV^^I < M^{P^x^) a.e. on fi} . (4.3) 

Set iy^'° = W^-^ G f/o^ £f = fX^ e S!" and qf = q^^^ G 5\ where 
we choose 0° = Q^ = 0. 

For m > 1, gh^en iterates W^'""^ G U^, ff"-^ G S!" and ^f"-^ G 5^ 
then 
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(i) Find W)^"^ e Ul; such that 

(4.4a) 
(ii) Find £/" e R''{W2n such that 

V V;'^ e g'iWl'"') . (4.4b) 

(iii) Set 

9.7"" = 9T~^ - P (VVr^'" - £1"") ■ (4.4c) 

Step (i) is equivalent to finding the unique W"^"^ G Uq solving the linear 
problem 



.v'] +p (viy;''" - fr .^vi - (Qr ' ^^ 

(r,^^') VT^'ef/o". (4.5) 






Step (ii) decouples to solving the problem on each element a E T^. Let 

^pla : f/o'' U x^'^ X ^'^ ^ ^ be such that 

Ci'-iv^t^^') = E ^i:M,£^v!:.) ' (4-6) 

where the subscript a denotes restriction to the element a. Hence, for all 
a eT , first find d) , G M such that 

—A, a 

4;:(WAr, iZ^ qTaT') ^ <'-"(^Ar' ^' qTaT^ v a g M^ (4.?) 

then project ' to the ball of radius [M^ {P'^W2'"^)]a- centred at the origin 
to yield 0"'*". The minimization (14. 7p leads to 

^ _ _^A^ . ^^_g^^ 



-A,(T 
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and we then set 



-A,a 




{M^(P''W2-'")], 


otherwise. 







< [M^{P^W'^"')], 



(4.8b) 



So overall, f/^ G g^W^"^) is such that 



{p{^ 



n,m 

iLa 



^WY^) + Q"''"-\ ^"^ - 0"''") > V V^'^ G ^'^(Vr^'™) . (4.9) 



On noting that 0" = VWj^, n = 1 — )• A^, in the variational case, Me(-) = 
fco £ IR>0) then following the abstract framework in §5.1 in [1^ one can show 
that for n = 1 — )■ A^ and m> \ that 



\Q\ 



'n,m— l,o o ,~-?i,m— l,r, 



> IW"'™''^ 



2 lv7T?/".m '7"'"' 1 2 






-71, m, 



2 

A lO.f^ 



2 I Z |V-7TTA","1 ^ ' |Z I 2 I i ' 

A \h^P IVW^A -^A lo,n+P \<tA 



'n,m—l , 



J-A 



lo,n 1 



(4.10) 



where PF^'™ 



w^A - wr, £r 



II - £^- and q7 



q:;-q; 



Hence, one can deduce from (14.101) . (I4.4cp . f l43|) and f lOj) for n = 1 -> A^ 
that as m — 7- oo 



^n,m ^ p^n^ 



f^- -> 0;; 



viy^. 



^r ^ ^A ' 



(4.11) 



Although we have no convergence proof of the iterative algorithm (I4.4al -c) 
in the quasi- variational inequality case, in practice it worked well. 



4.2 Solution of (Q^'^) 



Adopting the notation (I2.56p . we find Q" solving (I2.57p . and hence {W^^., 

Q^ } solving the n^^ step of (Q^^), (I2.54al b). using the following iteration: 

Set Q'f = Ql'^ G V^. For m > 1, given iterate Q"'"*"^ G V'', find 

gnm ^ yh ^^^-^ ^-^^^ 
-^B,T — 

(M,'^(^"-r„V.g 



I^R.. 1(5 V 






■5 -i-_B,r * 



rn(Y.Qr'V.t; 



Ii-_B,r 
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l'h/rh/,,n ^ v7 /^"l"l^l^ \/^n,m—l\r—2 \ /^ri,m—l\T—2 /^n,m—l „.h\h 



'-B,r '0 '2iLB,r 



'-B,r 



+ (^",V.t;^) \/v^eV_\ (4.12) 

where l^^l^ := (|wp + (5^)2 with 5"^ <^ 1. Clearly, the linear system ( I4.12p 
is well-posed. We note that similar algorithms have been used in [5l [6]. 
Although we have no convergence proof of fl4.12p . in practice it worked well. 

5 Numerical Experiments 

In this section we perform numerical experiments for our finite element ap- 
proximations (Q^'"^), f l2.22at b). and (Q^'^)) (I2.54at b). as stated in Section m 
except for ease of implementation we replaced Wq and Wq , fl2.15p . by wq 
and Wq = P'^[tt^Wo], respectively, in M^(-), f l2.16p . and in the initial data for 
both approximations. 

The approximation (Q^'^) is simpler and is easier to implement than 
(Qb^)' which is based on the lowest order Raviart-Thomas element. We refer 
to [3] for a Matlab implementation of the lowest Raviart-Thomas element. 
Although, both approximations lead to an efficient numerical approximation 
of the evolving sand surface, our numerical experiments, see below, show that 
only the (Q^'^) approximation yields a useful approximation to the surface 
sand flux for a reasonable choice of discretization parameters. 

The simulations have been performed in Matlab R2011a (64 bit) on a PC 
with Intel Core 15-2400 3.10Hz processor with 4Gb RAM. 

In all of our examples, we set the sand internal friction coefficient kg = 0.4 
and chose r = 1 + 10"'^ for the approximation (Q^'^), 02.54a[ b). 

The stopping criterion for the splitting iterative algorithm, (I4.4al -c). for 
(Q^'^) was chosen as 

\\TjTn,m Ti^n,m— 111 \\ in,m in,m—l\\ 

(5.1) 



, A \\L\n) 110; ll[Li(f7)]2 



For solving the nonlinear algebraic system, arising from (Q^'^) at each 
time level, we chose 5 = 10^^ for the iterative method (14.121) with stopping 
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criterion 



Ei^ii^ 



B,r (e) - Qb. (e) 



ee^'" , o 1 r^-4 






<3.10-\ (5.2) 



B,r vej 



Here £^'^ is the collection of edges associated with the partitioning T^ so that 
any v^ G V^ can be written as v^{x) = ^e^£hV^{e) (j) (x), where {0 jeef' 
are the standard lowest order Raviart-Thomas basis functions, see [3]. 

5.1 Approximation (Q^' ) 

We start with a simple variational inequality example. Let sand be disposed 
onto a flat, Wq = 0, open circular platform, Q = {x : \x\ < 1}. The source / 
is uniform in its support \x\ < Rq = 0.2 with J^ f{x, t) dx = 1 for all t > 0. 
Due to the radial symmetry, the analytical solution to this problem, {w,q}, 
is easy to find. The growing pile starts as a cut-off cone having critical slopes, 
volume t and height t/(7ri?g). Then, at t* = vr A;o i?o a/3 ^ 0.0174, the pile 
turns into a cone w{x,t) = ko max(_Rc(i) — |ai|,0). This cone grows until its 
base, a circle of radius Rdt) = {St/{n ko))^, fills the domain Q. The flux can 
be found as q{x,t) = q{\x\,t) x/\x\, where q{R,t) is a solution to the balance 
equation 

;^^(^g) = /-^ fori?e(0,l), g(0,t)=0, fort>0. (5.3) 

The iterations of the augmented Lagrangian method with splitting, recall 
Section 14.11 converged quickly with p = 1. For t = 0.1 we compared our 
numerical approximations obtained for different finite element meshes and a 
constant time step r with the analytical solution. The approximate surfaces, 
W2 with nr = t, were close to the exact surface, w{-,t); see Figured! left. 
We checked that, for the meshes employed and r G (0,0.01], the error in 
w was dominated by the spatial discretization. For meshes with maximal 
element sizes h =0.01, 0.02, 0.04 the relative errors oiw{-,t) in the L^ norm 
were, respectively, 0.3%, 0.9%, and 1.8%. 

Although in our simulations the approximate flux iterates Q^"^ also con- 
verged on every mesh, no pointwise convergence of Q^, with nr = t, to the 
analytical flux g(-,t) was observed; see Figured! middle. The approximate 
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flux Q^ has a fine structure in the region where the exact flux is not zero. 
There elements a G T^ with zero numerical flux were intermixed with ele- 
ments in which the numerical flux was much stronger than the exact one; 
see Figure [H middle and right. Such a behavior of the numerical solution 
does not contradict our proof of its vague convergence to g; but, clearly, a 
different method should be employed for approximating the flux even in the 
variational inequahty case. 
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Figure 1: Variational inequality, {Q^) approximation, simulation results 
for h = 0.04, T = 0.01, t = 0.1. Left - exact surface w{\x\,t) (red line) 
and its approximation W^ at the mesh nodes (black dots). Middle - exact 
flux modulus li^da;!, t)| (red line) and its approximation \Q'^.\ in the elements 
(black dots). Right - levels of |(5"| showing the mosaic structure of the 



approximate flux Q" 



lliLA 



The situation is similar for the quasi- variational inequality case, with only 
the evolving pile surface being approximated well using this method. In our 



second example, see Figure |21 we set wq 
— (0.3, 0), for the square fi = (—1, 1) x 



max(0.5 — 



X 



X 



.Oh 



0), where 



(—1, 1). The source / is uniform 
1 for all t > 0. For h = 0.02 the 



in its support |x| < 0.7 with J^f{x,t) 
generated mesh contained approximately 34,000 elements. Since the gradient 
constraint is now updated after each iteration of the splitting algorithm, 
existing theory does not guarantee its convergence. We found that, for the 
regularization parameter e = 0.01, good convergence of this algorithm is 
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41^ <#■• #^ 



Figure 2: Regularized quasi-variational inequality with e = 0.01, (Q/^) ap- 
proximation, simulation results for h = 0.02, r = 0.01. Left - the support 
surface Wq, the dashed line indicates the boundary of the support of /. Mid- 
dle and Right - approximate sandpile surface, W^, corresponding to t = 0.1 
and t = 0.2, respectively. 

achieved for a smaller value of the augmented Lagrangian parameter p. In 
this example we chose p = 0.05 with stopping criterion (15. ip and obtained 
the solution with r = 0.005 and twenty time steps in 11 minutes of CPU 
time. 

5.2 Approximation (Q^'^) 

This approximation performed better in the variational inequality example 
from the previous section, recall Figure [1] Using the time step, r = 0.005, 
and meshes, h=0.02 and 0.04, we obtained, for t = 0.1, the pile surface 
with smaller relative errors in the L^ norm; 0.1% and 0.6%, respectively. 
Furthermore, for this approximation the fluxes Q'^ ^ also converged to the 
exact solution. Comparing the approximate and exact fluxes at the element 
centers we estimated the relative flux error in the L^ norm. For the two 
meshes chosen these errors were, correspondingly, 2.8% and 5.2%. 

Choosing a different initial support, Wo{x) = max(0.4 — |a;|, 0), and keep- 
ing the same source / and domain Q from the variational inequality example, 
we arrive at a quasi-variational inequality problem that can be solved an- 
alytically for the unregularized M(-). Being discharged from the source, 
sand now pours down the steep conical part of the support surface and 
forms a pile around this cone. The volume of the pile is t and its sur- 
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face w{\x\,t) = kQ{R2(t) — \x\) for \x\ G [Ri(t),R2(t)], see Figure El Using 
simple geometric arguments, we first find the two variables, Ri{t) and R2{t), 
determining this surface from the equations 



t = ^ [{Rl - Rl) ko - (0.4^ - Rl)] and R2 



Ri + —{OA-Ri] 



We then find the flux using the balance equation ( 15. 3p . 

We solved the problem numerically, see Figure HJ with the regularization 
parameter e = 0.005 and estimated the errors of WJ^,, and Q" at t = 0.1 
using the analytical solution. As could be expected, the smaller the value 
of e, the more difficult it is to obtain convergence of the iterations f l4.12p 
in the quasi-variational inequality case. We were, however, able to achieve 
convergence of these iterations by decreasing the time step r. For the stated 
value of e, we chose r = 0.0005 yielding 200 time steps on the time interval 
[0,0.1] for two different meshes. For a mesh generated with h = 0.04 the 
relative errors in the L^ norm were 0.6% for the pile surface and 5% for the 
surface flux. For a finer mesh, h = 0.02, the corresponding errors were 0.1% 
and 2%. These results confirm the validity of our regularization, M^{-), of 
M(-). 




Figure 3: Sandpile, yellow region, forming upon the support platform with a 
steep cone. 



^/l,r ^ 



We solved again, now using {Q^\), the quasi-variational problem consid- 
ered in Figure [2] above, using the same mesh, time step, and the value of 
regularization parameter. Now we were able to find good approximations. 
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a 0.2 





Figure 4: Regularized quasi-variational inequality with e = 0.005, (Q^'^) 
approximation, simulation results for h = 0.04, r = 0.0005 and t = 0.1. 
Left - exact surface w{\x\,t) (red line) and its approximation W]^^ in the 
elements (black dots). Middle - exact flux modulus |g(|a;|,t)| (red line) and 
its approximation \Q^ \ at the element centers (black dots). Right - the 
Qo vector field, where the dashed line indicates the exact pile boundary. 



not only to the pile surface but the surface flux as well, see Figure [5l and 
the computation time was about the same. We note that as the surface H^^r 
touches the support boundary dfl at some time in the interval (0.1, 0.2) sand 
flows out of the system, which can be seen from the flux Q^^. 

In our last example Q = (—1,1) x (—1,1) and wo = min(max(|a;i| — 
0.9, |x2| — 0.9), 0) is the surface of an inverted pyramid supplemented, to 
satisfy the no- influx condition (II. 8p . by a narrow horizontal margin. The 
uniform source is f{x,t) = 0.25. Sand, discharged from the source, flows 
down the pyramid faces until it reaches a pyramid edge; then it pours down 
along the edge and forms a pile above the apex of the inverted pyramid, see 
Figure El Our numerical solution clearly shows the singularity of the edge 
fluxes. 



48 















Figure 5: Regularized quasi-variational inequality with e - 
[21 (Q^^) approximation, simulation results for h = 0.02, r 



0.01 as in Figure 
: 0.01 and t = 0.2. 



Left - the calculated surface WJ^^. Middle - the flux modulus |Q" | at the 



element centers. Right - the Q" vector field and level contours of Wj 



-il_B,r 



B,r- 






Figure 6: Regularized quasi-variational inequality with e = 0.02, (Q^^) ap- 
proximation, simulation results for h = 0.02, r = 0.0025 and t = 0.075. 
Left - initial surface wq (blue lines) and the approximate surface VTg ^ in the 
elements (grey surface). Middle - the flux modulus \QZ I at the element 



centers. Right 



91, 



lO" I 



vector field and levels of W^ , 
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